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Abstract. The ATLAS detector is capable of resolving the highest energy pp collisions at lumi- 
nosities sufficient to yield lO's of simultaneous interactions within a bunch collision lasting <0.5 
nsec. Already in 2011 a mean occupancy of 20 is often found in pp running. In 2004 studies by 
ATLAS showed that the detector would have excellent performance also for the foreseeable parti- 
cle multiplicities in the highest energy p-Pb and Pb-Pb collisions that the LHC will produce. These 
studies resulted in a letter of intent to the LHC committee by ATLAS to do physics with these beams 
also. In the past 2 years of data taking, ATLAS detector performance studies have confirmed these 
expectations at the actual multiplicities presented below. 

The ATLAS program removes an artificial specialization that arose about 30 years ago in high 
energy physics when the energy and intensity frontier moved to colliders. Before that time, for 
example, the same experiment that discovered the T (CFS and E605 at Fernulab) also measured the 
nuclear modification factor in the production of high pr identified charged hadrons using nuclear 
targets from Beryllium through Tungsten. 
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INTRODUCTION 

In roughly 1974 interest began to develop in studying the highest energy hadron-nucleus 
and nucleus-nucleus collisions. That year is sometimes identified w^ith a triumph of e^e~ 
colliders- the "November Revolution" [1] as well as the first papers on QCD[2]. Also in 
the same year a paper by Lee and Wick asked "w^hether it is experimentally possible 
in a limited domain in space to excite the ordinary vacuum to an abnormal state" [3]. 
A conference in Bear Mountain, NY, organized by Lee, Lederman, Ruderman et al.[4] 
explored the experimental possibilities, including using an existing collider (ISR) with 
nuclear beams. In the same year interest developed in hadron-nucleus collisions as a tool 
to understand the space-time development of particle production due to the Fermilab 
experiments of Busza[5] et al. Since the LHC provides access to both the highest energy 
pp collisions and the highest energy collisions of nuclear beams it is natural to use the 
largest and most complex particle detectors (e.g. ATLAS) to explore these topics also. 

ATLAS can also bring about a convergence with the proton and nuclear structure 
studies traditionally done in electron-nucleus (or proton) collisions (i.e. HERA). Both 
the p-Pb data[6] and the hard photoproduction data[7] will provide unprecedented access 
to low-x structure of hadrons. 
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FIGURE 1. The x and y components of missing Et resolution vs. Et showing that the same behavior 
holds for Pb-Pb data as in pp. 

ATLAS DETECTOR AND DATASET 

The primary focus in preparing ATLAS for Pb running was to ensure good track re- 
construction efficiency and small fake fraction for charged tracks in the inner detector 
for the highest multiplicity events. This dominates the track quality for measurement 
of charged multiplicity, nuclear modification of the charged particle pr spectra as well 
as for the reconstruction of muons (since the outer tracking region for muons always 
has low occupancy, as expected). The inner tracking efficiency was well modeled by 
simulations resulting in an uncertainty of 3%. 

The ATLAS calorimeter is highly segmented, both longitudinally and transversely 
and covers the full rapidity region to ±4.9. Nevertheless in the highest multiplicity 
Pb-Pb events the underlying event, due the high multiplicity of soft particles, will 
typically contribute to the energy deposit in a given tower being used to measure jet 
energy, for example. The underlying event subtraction must be determined from the data 
themselves, as discussed below. It is important that the calorimeter response vs. energy 
deposit be well understood. This is demonstrated in Fig. 1 , where the x and y components 
of the missing transverse energy resolution are plotted for both pp collisions (which have 
lower transverse energy than is possible in Pb-Pb) and for Pb-Pb collisions. 

In the 2010 Pb-Pb run ATLAS recorded approximately 9 iib~^ of luminosity at a 
recording efficiency of ^ 95%. Of these data, roughly 10% were recorded with central 
magnetic field turned off in order to increase the efficiency of low momentum tracking 
for the charged multiplicity measurement , discussed below. The minimum bias trigger 
consisted of the Minimum Bias Trigger Scintillators ( MBTS, covering the region 
2.09 < |r;| < 3.84) since these were familiar from the pp runs and, for most of the 
data. Zero Degree Calorimeters which subtend an angle of ^ 350 /irad with respect to 
the beam direction and accept beam energy neutrons from nuclear breakup out to pr ^ I 
GeV/c. 
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FIGURE 2. Charged multiplicity density per participant nucleon pair for the most central ATLAS 
collisions (top). Dependence on the number of participants (bottom). 



All data were collected with a beam energy of 1.38 TeV/nucleon and a typical lumi- 
nosity of a few X 10^^ yielding an inelastic collision rate of lOO's of Hz. 

MEASUREMENT OF CHARGED PARTICLE MULTIPLICITY 

AND CORRELATIONS 

Charged particle multiplicity and correlations are more difficult to interpret in pp data 
than in Pb-Pb because, up to now, the experimental tools for characterizing the initial 
conditions of pp collisions have been limited. In contrast, for the Pb-Pb data presented 
here the initial collision impact parameter (or, by calculation, the number of participant 
nucleons ) and orientation of the collision are determined with good precision. 

In the case of pp data, for example, multiplicity is measured for a sample of collisions 
ranging from peripheral (and, in some cases, diffractive) to very central for which we 
would expect more or less copious particle production. Similarly the leading baryon 
effect implies that a large fraction of the collision energy is carried away by a leading 
baryon (protons or neutrons with roughly equal probability) in each hemisphere. It has 
been shown that when the effective c.m.s. energy of pp collisions is corrected for the 
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FIGURE 3. Vn vs pT for several centrality classes. 



leading particle effect, pp multiplicities exhibit the same trend vs. energy as those from 
e^e~ or even nucleus-nucleus [9]. 

In the case of Pb-Pb collisions the impact parameter is determined by mapping the 
frequency distribution in any sensitive variable into the known frequency distribution of 
collisions vs. impact parameter, ie: 

Luminosity (b) oc bxdb (1) 

In the present analysis the frequency observed in ZEt in the forward calorimeter is used. 
The charged multiplicity is then renormalized by the calculated number of participant 
nucleons for this collision [10]. 



CORRELATIONS 

In Pb-Pb collisions the azimuthal distribution in either the forward or backward (or both) 
hemispheres can be used to learn about the initial geometry of the collision. Because the 
reaction plane (i.e. the orientation of the impact parameter vector, b) is independently 
measured in each hemisphere, the resolution of the reaction plane can be determined 
directly from our data. One method to measure the reaction plane (full FCAL method) 



consists of determining the plane in each forward calorimeter (3.2< |r; | < 4.9) from the 
energy distribution in the calorimeter cells of the first 2 depth layers and combining the 
2 measurements. The azimuthal distribution of charged particles measured in the inner 
tracker can be expressed as a Fourier series, ie 

^ oc LvnCos{n^) (2) 

where O is the azimuthal angle of tracks relative to the reaction plane. The first 
moment of the distribution, called directed flow must by symmetry, equal zero at central 
rapidities so it is not presented in the analysis here from the inner tracker. We find that 
the flow moments vary little with pseudorapidity but depend strongly on track pr and on 
centrality (given by equation 1). The panels in Fig. 3 are divided into centrality classes 
with the upper left taken from the 5% smallest impact parameter events [1 1]. 

The flow moments, v^, are sensitive to the hydrodynamic properties of the hot, dense 
precursor to hadron production in these events. One hydrodynamic calculation [12] at- 
tempts to account for the relative strength of the Vn for various choices of viscosity of 
this medium and obtains approximate agreement. 
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FIGURE 4. jet yield for various centrality classes divided by yield in a peripheral class. R=02, 0.4. 

Hard scattering of partons within Pb-Pb collisions can be used to probe the field of 
high parton density within which the scattering occurs since the hard scattering occurs at 



the very beginning and the scattered parton must traverse the medium before it emerges 
and fragments into a jet. There are a number of ways that the parton could lose energy 
as it traverses the medium, including radiative and coUisional energy loss. Although 
fragmentation occurs late in time compared to the time it takes the parton to traverse this 
hot medium, gluons radiated during the energy loss process would effectively modify 
the particle distribution (fragmentation functions) of particles in the reconstructed jet. 

If the energy loss is purely the result of elastic collisions the jet fragmentation could 
look like that of a jet produced in a pp collision but the reconstructed jet axis might be 
deflected. This would lead to a modified distribution in the ^jet-Jet correlation of back- 
to-back jets. In fact, ATLAS sees small modification of the back-to-back correlation of 
di-jets in central collisions (shown elsewhere [14]). Studies of jet modification in central 
Pb-Pb collisions could shed light on the process of vacuum fragmentation, which is 
essentially treated as an empirical parametrization in most of high energy physics. 

Jet reconstruction in heavy ion collisions is more difficult in Pb-Pb data than in pp and 
doing it correctly has been a focus ever since the original studies for the ATLAS Heavy 
Ion LOI. As has been pointed out in the literature, incorrect treatment of the underlying 
event energy and its fluctuations can easily lead to artifacts [13]. 

The method used by ATLAS to reconstruct jets in this environment attempts to do 
this making the best use of the detailed measurement of energy flow in ATLAS' large 
acceptance and highly segmented calorimeter. 
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FIGURE 5. Longitudinal (left) and transverse fragmentation distributions for various centrality classes 
divided by a peripheral class. R=0.4 

For results presented here jets are reconstructed using an mti-kf algorithm and cone 
radii of both R=0.4 and R=0.2. Layers of the calorimeter within a tower of 0.1 x 0.1 are 
combined taking into account energy-density dependent corrections arising from non- 
compensation of the calorimeter. We then calculate a distribution of energy flow from the 
underlying event, based on energy deposit excluding the jet window. The modulation, 
due to flow of particles in the underlying event is determined using the results presented 
above and the reaction plane determined with the full FCAL method. This method, 
nevertheless, slightly overcompensates since a fraction of the jet energy contributes to 



the normalization of the underlying event energy. A second iteration of the subtraction 
procedure is therefore performed. Results presented here are for jets with \rijet \ < 2.8. 

Fluctuations of the underlying event (tower-by-tower and for different layers) are 
obtained from simulation with HIJING+GEANT and compared to data. Our conclusion 
is that fluctuations in the data are reproduced to within 10%. 

Complementing an earlier publication which showed evidence of jet quenching using 
di-Jet events we have measured inclusive jet production and compared the yield in sev- 
eral centrality bins to the yield in peripheral (80%) collisions. Figure 4 shows the yield 
expressed as a ratio of central to peripheral. The suppression reaches ^ 0.5 independent 
of jet energy for both R=0.2 and 0.4 jet cone radii. This is a similar suppression to what 
was observed for high pr charged particles at RHIC, which leads one to suspect that the 
fragmentation function for jets produced in central Pb-Pb collisions is unmodified from 
the vacuum distribution. 

Fig. 5 shows that this is true in detail for the case of R=0.4 reconstructed jets. Similar 
results are obtained for R=0.2[14]. 

J/1// SUPPRESSION AND REFERENCE MEASUREMENT 

Since the jet yield is modified so dramatically it would not be surprising to see suppres- 
sion of a di-quark bound state, such as the J/i// . In fact this has been observed early in 
the SPS Heavy ion program by NA50. There are many aspects of the observed suppres- 
sion which are not well understood, however. Originally it was felt that this suppression 
must occur because color screening in Heavy Ion collisions would prevent formation 
of a quark-antiquark bound state. One way to get further insight into this problem is 
to search for similar suppression of or bosons since these are neither strongly 
interacting or composite. Any modification of their yield would have to arise from the 
initial state (i.e. due to nuclear structure function). 

ATLAS has both measured J/i/A suppression and searched for similar suppression of 
the weak vector bosons. 

Dimuons were selected with opposite charge, pr >3 and |t]| <2.5. After sideband 
subtraction to remove dimuons from the continuum a sample containing 190± 20 J/i// 
candidates was obtained for the most central 10%. A suppression in this most central bin 
relative to a peripheral sample of ^ 0.5, is found which is not very different from the jet 
suppression result. 

To obtain a sample of candidates muons were selected by removing events with 
a high mass dimuon (m^^ > 66 Gev) and the inclusive spectrum was fit to combination 
of the continuum spectrum -mostly due to charm and b quark decays- and a W-decay 
template tuned on a W sample from pp data. The inclusive spectrum was then fitted to 
obtain the number of W candidates in a given centrality class of Pb-Pb data, as illustrated 
in Fig. 6. Also shown in Fig. 6 is the yield ratio for peripheral collisions divided by the 
yield in a given centrality class, which shows consistency with what you would expect 
from binary collisions scaling- i.e. no suppression. 
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FIGURE 6. (left)Inclusive muon pr distribution with cuts described in the text. Also shown is the fit 
to W decay template and the background from heavy flavor decays, (right) The W yield vs. centrality 
selection presented as peripheral yield divided by yield in the centrality bin, corrected for the number of 
binary collisions. 



OUTLOOK 

In November 201 1 the LHC will again collide Pb beams (yielding an expected ^ 4-fold 
increase in integrated luminosity) and for a short test, possibly yielding of order 10^ 
inelastic collisions, will attempt p-Pb collisions also. 

There are a number of further measurements that could be used to clarify the mecha- 
nism of jet quenching. CMS has shown that the energy loss that results in jet quenching 
appears in the form of low pr particles at wide angles to the jet. The anisotropic emission 
of particles described by the Vn moments is modeled as small initial state perturbations 
which propagate hydrodynamically. It would be interesting to study the Vn for a class 
of events where a ^ 50 GeV jet is emitted since its energy loss must produce a large 
perturbation. It would also be of interest to compare jet suppression for b-quark jets. 

Near r;=0, d-Au data at RHIC showed an enhancement rather than a suppression of 
high pt charged particles as did E605[8]. However at forward rapidities a suppression 
was observed. New data on p-Pb will be useful for discriminating among present mod- 
els for this suppression. CMS has shown that similar features of correlations in Pb-Pb 
data also appear in pp collisions selected for high multiplicity which likely corresponds 
to very large gluon densities[15] and studies are under way which reveal other similar- 
ities to high multiplicity Pb-Pb collisions [16]. It would be interesting to explore other 
methods [17] to select pp collisions with similar properties. 
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